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Abstract 

On| , In [1] a linear method of solving a particular set of Lichnerowicz-type equations through 

the implicit function theorem was sketched in order to implicitly construct Shape Dynamics' 
global Hamiltonian and eliminate second class constraints. This method was completely laid 
out in [2], and in [3] it was used for extending Shape Dynamics (SD) to the non- vacuum case, 
showing how other fields are coupled to the theory. In the latter paper it was noticed that, 
unlike the vacuum case, the use of such methods yielded puzzling bounds on the density of 
some types of fields. Here we show that the original SD cannot be extended beyond such 
■ bounds, but that a slight modification of the original can withstand any type of coupling. 

When the bound is broken, the theory does not come equipped with a single Hamiltonian as 
in vacuum SD, but with a finite set of weakly commuting Hamiltonians, which we describe. 
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^Jj 1 Introduction 



^ | 1.1 Shape Dynamics in a nutshell 



Let us start by very briefly defining the setting in which this study takes place: the theory of 
\Q . Shape Dynamics (for a more comprehensive introduction see [4, 1, 2]). 1 In very brief terms, 

0^ | SD is a classical theory equivalent to pure gravity in which the spacetime picture that underlies 

GR is replaced by a picture of evolving 3-geometries, as in the original intention of Wheeler. 
What makes SD unique amongst more standard geometrodynamical theories of gravity, is that 
| the usual invariance under spacetime refoliations is traded for invariance of the theory under 

local spatial conformal transformations that preserve the total spatial volume. The equivalence 
is obtained by manipulating the constraint structure of General relativity, but can only be 
achieved if one first uses the Stuckelberg mechanism to extend the original phase space with 
5_i ' extra degrees of freedom. After this "Stuckelberg extension" of the original ADM dynamical 

system, one obtains a system even more redundant in symmetries, 2 which is reducible to both 
the original ADM and SD through the use of different gauge conditions on the added variables. 

The basic condition required for the symmetry trading is the existence of two (the original 
and the "new") symmetries which are maximally symplectic (i.e. gauge fix one another). Al- 
though this basic condition can be expressed in the original ADM phase space, the construction 
of SD itself requires us to use extended phase space explicitly. It is then the implicit function 
theorem which enables the use of a linear method and the maximally symplectic character of 
the dual symmetries to obtain SD from extended phase space. This linear method, fully laid 
out in section 4.3 of [2], works perfectly well in vacuum SD. It shows that any ADM vacuum 
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For the original motivation to regard the role of conformal transformations in ADM, see [5]. 
2 In fact, the most technically clear way of obtaining SD requires first of all a trivial embedding of the original 
system into the extended phase space, and only then a canonical transformation that in effect introduces a 
conformal character to the system. 
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dynamical system has its SD counterpart, provided we interpret SD in a purely geometrody- 
namical way, as explained here in section 2.2.2. However, as the author found out, the original 
linear method was not comprehensive enough to deal with certain non- vacuum systems. This 
result was first reported in [1]. It is with the extension of this method that the present paper 
is concerned. 

1.2 Non-uniqueness of the SD Hamiltonian 

The original SD construction allowed the linear method to reduce an extended theory to a theory 
existing on the ADM phase space, possessing a single global Hamiltonian and local conformal 
invariance. It relied heavily however on the positivity of the linear term in a differential operator 
arising upon reduction of the system. Unfortunately, once one tries to extend the result to couple 
matter, the required term is no longer positive, and thus the relevant differential operator may 
have non-trivial homogeneous solutions. Conceivably, these could have yielded different theories 
of Shape Dynamics, or perhaps the reduced system could have turned out to be inconsistent. 

The present paper will use ellipticity of the relevant differential operators to extend that 
result and show that in the general case one still retains local conformal invariance, but must 
allow for the possibility of a finite set of commuting global Hamiltonians instead of the single 
one. The number of commuting global Hamiltonians is related to the kernel of the relevant 
differential operator. Here we show that there is still a coherent "Shape Dynamics" theory 
arising, at least for the case of a compact closed 3 manifold. 

Before we start, we give a warning about the level of mathematical rigor aimed for in this 
paper. Although the paper contains proofs and propositions, it is aimed at a physics audience. 
Its purpose is to arrive at a new physical result regarding SD, namely that for some types of 
matter (including certain values of the cosmological constant) we have a theory dynamically 
equivalent to ADM with a finite number of compatible notions of "absolute" (global) times. We 
are not extremely concerned with the type of functional spaces we will be dealing with, and 
often sloppily characterize distributions as elements of C°°(M), and other such mathematical 
atrocities. However, we have aimed to give robust, plausible arguments, such that if due care was 
given to such domains all statements here could be put in a firm functional analytic grounding. 3 
We will usually assume that we are dealing with Hilbert spaces endowed with an L2 inner 
product. 

2 The vacuum case 

We now give a more technical, but very streamlined account of the construction of SD. This will 
be useful to understand the limitations imposed by the bound, and the attempt to overcome 
these later. 

2.1 Technical setting of SD 

Let us now briefly review the setting for the construction of Shape Dynamics as a theory 
equivalent to ADM gravity on a compact Cauchy surface S without boundary. We start with 
the standard ADM phase space r ADM = {(g, tt) : g G Riem, it G (Riem)}, where Riem denotes 
the set of Riemannian metrics on the above defined 3-manifold E, and the usual first class ADM 
constraints, i.e. the scalar constraints 

n ab _ 12 

S(x) = ^ 2 " - V~9R (1) 

3 For example one could use Sobolev completions to given finite rank Sobolev spaces and then inverse limits 
to attain such rigour. Or simply work with C 2 (M) instead of C°°(M). 
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and momentum constraints 

H a (x) = n%{x) (2) 

thereon. We extend the ADM phase space with the phase space of a scalar field 4>{x) and its 
canonically conjugate momentum density 7T^(x), which we introduce as additional first class 
constraints Q(x) = 7r^(x) ~ 0. The system is thus merely a trivial embedding of the original 
ADM onto the extended phase space r ADM x and no sign of a "conformal transformation" 
is in sight. 

It is by introducing the canonical transformation 7^ generated by the generating func- 
tional F = J d?x [g ab e^Ii ab + (pU^j , where <p(x) := <p(x) - \ ln(e 6( ^)g using the mean (f) g := 

y J d 3 X\JJg\f(x) and 3-volume V g := J d?x\J\g\ that a conformal character of the theory starts 
to appear. For it is this transformation that emulates a volume-preserving conformal transfor- 
mation in the original canonical variables: 

9ab{x) -> T^Qabix) : = e 4 ^ g ab (x) 

7r ab {x) -> T^7r ab (x) := e" 4 ^) (n ab (x) - ^^(vr)(l - e 6< ^)) ^ 

4>{x) ->■ T<f>(p(x) := 4>(x) 
■K(j>(x) ->■ T^ir^(x) := tt^x) - 4(tt(x) - (it) y/g). 

Renaming for convenience: 

4(tt(x) - <tt) ^5 =: £>(x) 
at this point we have the first class set of constraints 

T<j>S(x), T<j,H(x) and T<j>Q(x) = tt^x) - D(x) (4) 

We should note however that the extended system does not yet possess any notion of "conformal 
symmetry", for the action of T$Q{x) on the canonical variables (f) and tt^ is not conformal. We 
also note that 7^Q(x) acts trivially on phase space functionals of the form T,p(f(g, vr)(x)) and 
that T^Hix) acts (weakly) as diffeomorphisms in extended phase space. 

It is the presence of the volume-preserving element of the conformal transformations that 
accounts for the more complicated structure of the momentum transformation, and it is also 
the source for the presence of a non-trivial global Hamiltonian in SD, as we show now. 

To regain ADM, one imposes the gauge- fixing <f> = 0. To obtain SD, one introduces into 
the extended system above the gauge-fixing ir^ = 0. To get SD from the imposition of this 
gauge-fixing is where it is most convenient to use the linear method, as we now explain. 



2.2 Finding and solving the second class constraints 

Again, this is a streamlined version of the entire argument. For a more comprehensive and 
in-depth treatment, see [2], section 4.3, where the method was first completely laid out. As 
a matter of nomenclature, we say that a system of constraints is first class if it commutes on 
the constraint surface (weakly), and it will be "purely" second class, if the Poisson brackets 
are invertible (for example, if they are proportional to the identity). If this is the case the 
bracket itself cannot impose further restrictions on phase space (constraints), and the second 
class constraints must be either completely solved for (used as definitions of given canonical 
variables), or one must use the Dirac bracket to project the symplectic structure onto the 
constraint surface. We will show that the first case is attainable here in an elegant manner 
which gets rid of extra variables. 
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2.2.1 Separate out the first class component 

The only weakly non- vanishing Poisson-bracket of the gauge-fixing condition ttAx) = with 
the constraints of the linking theory is 



{T+S(N),*4(x)} = 4T<t,{S(N), tt(x) - (vr) Jg(x)}, (5) 

By the canonical transformation properties of the theory, we can translate any assertion made 
at <p = to a different <j) by the use of T^. Hence the relevant Poisson bracket becomes: 

4{S(N),tt(x) - <tt) ^g(x)} = 8AN(x)^ - 8(AN)^- 6S(x)N(x). (6) 

In (6) A is the differential operator: 

A = V 2 -^(n)n-R. (7) 

On the gauge fixed constraint surface T ( pS(x) = T$Q(x) = tt^ = at <j) = 0, 4 using the 
momentum split ir ab = o ab ^fg + \g ab ^ into its trace and traceless part, a ab = a ^/g, we have 

A^V 2 -^(7r) 2 -a ab a ab . (8) 

Thus our operator A can be written as A = V 2 — A where A[g, tt;x) := a ab a a f, + i (ir) 2 > 0, 
i.e. it is a positive-definite function, vanishing only when ir ab = 0. Rewriting (6) in a more 
convenient manner we obtain: 

-Us(A0, vr(x) - (vr) ^(x)} « 2(AiV - (AN)) (9) 

The first question to be asked now is: are the given constraints purely second class? I.e. 
is the Poisson bracket matrix (6) invertible? The answer is dependent on the existence of 
homogeneous solutions, N(x) ^ such that: 



AN 



(an) (10) 



If there is such a solution, the answer is no, the bracket is not invertible as it stands. 

The initial challenge then is to try to separate the constraints T^S into a purely second class 
part, which thus has invertible Poisson brackets with ir^ = 0, and is thus completely gauge-fixed 
by this condition, and another part, which weakly commutes with the gauge-fixing condition. 

From the positivity of the linear term A in A, we understand that it cannot have a non-trivial 
homogeneous solution on a compact manifold, i.e. AiV = implies N = 0. This follows from 
comparing both sides of AiV = at the maximum (minimum) of N. An equivalent statement 
is that for any function / there exists a unique Nf such that AiV = /. Furthermore, a solution 
of (10) has to be such that ANq = c, where c is a spatial constant. Since any constant will 
do, we have a linear space of solutions at each time, as it should be. We fix a generator of this 
space by adjusting the constant c to be such that (iVo) = 1 which is guaranteed to be possible, 
since for c < 0, N(x) > for all x. The truth of this last statement is apparent in the following 
manner: suppose that N(x) < for some x, then for some y, N(y) will attain a minimum, and 
thus V 2 N(y) > 0. The equation then becomes A = —S/ 2 N(y) — |c|, which is absurd. Note that 
almost all of the assertions we have made so far require positivity of A, and it is exactly this 
positivity that will be broken in section 3. 



4 We have taken all the relevant brackets, so now it is legal to regard the behavior of the operator on the 
constraint surface, and at (f> = 0, and then generalize by using the canonical transformation properties (which 
preserve the brackets). 
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2.2.2 Using the linear method to solve the second class parts. 

We have succeeded in separating out one generator of the first class part of 7^5 with respect 
to 7r^, namely T(j,(S(No)), where Nq obeys (10) and (Nq) = 1. We have not yet singled out the 
uniquely second class part, which has invertible Poisson bracket wrt tt^. After we have done 
that, to arrive at Shape Dynamics, we still must show that that second class part is solvable 
by a definition of the conjugate variable to tt^, i.e. that we can choose (f> = cp (g,ir) such that 
this solves the second class component of T^S. We will then have the complete theory on the 
original phase space of ADM, with constraints T^, (S(Nq)) , D , H. It is here that we will use 
the implicit function theorem. 

First, it will prove very useful to be more careful about what we mean with "a smearing". 
Consider 

TS(x) : T x r*(C°°(M)/V) -> C°°(M), (11) 

where V is the linear space of non-zero constant functions and the overline means we take the 
de-densitized version. I.e., suppose we have some density: F : T — > D(M), where D(M) is the 
space of densities of weight one over M. Then F(g,ir) = F(g,ir)/y/g. 5 

Since the map (11) does not depend on ir^, we can fix n^x) = f(x). Then 

TS(x\, =f[x) : r x C°°(M)/V -> C°°(M). (12) 

We note that in fact T<f>S(x) depends solely on C°°(M)/V. In any case, a smearing is given once 
we have established the L2 inner product: 

(F, h) := [ Fh^g~d 3 x = [ Fhd 3 x =: (F, h) (13) 
7s 7s 

This creates a pairing between the space of densities and that of functions, which we defined in 
the above rhs. 

Consider the linear operator: 

^7^(50,^0,4)^0 : T (C°°(M)/V) -> C°°(M) (14) 

Here we have denoted the derivative in the second coordinate, the one parametrized by 4>, by a 
subscript C, and where T X M denotes the tangent space to M at x G M, and, as in usual partial 
derivatives, one holds the coordinates (g,TT,Tr^) fixed. We will omit from now on the "initial" 
point (go, ttq, 7r°) where we take the derivative. It can be checked that the space Tq(C°° (M) /V) 
can be redundantly parametrized by functions / G C°°(M) as f —(f), so that constant functions 
are mapped to the zero element in Tq(C°° (M) /V) . 
The tangent map is given by: 



6 cT^S {(j)=0 := s ^ = {7^5(x),7r^(i/)}|^= = y/g(x)(A(x)5(x, y) — ) (15) 



Contraction of (15) with N(x) yields AiV — (AN). We note that contraction in the x variable 
requires us to use the adjoint 5cTif,S*^_ Q , and as this is not a self-adjoint operator the distinction 
is important. Thus 

(S C %S)* ■ N = AN - (AN) (16) 
and from the uniqueness of the kernel of the adjoint we in fact have 



Im(5cT^S)(BN = C oo (M) (17) 



5 We again emphasize that these statements could be made rigorous using the constraints to be taken as test 
functions and the smearings as the space of linear functionals over C°°(M). This would resolve some issues to 
come as for example the use of the delta function as a smearing. 
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Now we construct a modification of S(x), which we will call L such that it has the same 
tangent map but its range must be such that: 

(W^=/(x) : r x C°°(M)/V -> lm(6 c T*S). (18) 

This map is given by L = S(x) — S(No)^/g, since as one can readily check the tangent map 
Sc%f,L indeed stays (weakly) the same and {L,Nq) = (since (No) = 1). 

We have not yet proven that ScT^L is a topological linear isomorphism (or alternatively 
that TtfiL is purely second class with respect to ir^). We have shown that it is a surjective 
linear map (or alternatively that Kev(5cT ( j > L)* = 0), but we must still prove injectivity, i.e. that 
lm(5cT<f>L)* = Tq(C°°(M)/V). As we know, the elements of the latter space are given by / — (/) 
where / E C°°(M). The differential operator A is invertible, possessing a Green's function. 
Thus for any function / there exists a unique Nf for which ANf = f. Thus (dcT^S)* ■ Nf = 
ANf — (ANf) = f — (f), and we have surjectivity of the formal adjoint. Let us here briefly show 
an alternative proof of injectivity, which will be more useful in the generalization of the following 
section. First, we have that since we are assuming a Hilbert space structure, Ker(5cT(j>L) = 
(Im(ScT(/)L)*) . Thus, by the invertibility properties of A, any element of (lm(5cT ( pL)*) 1 - must 
satisfy (/ — (/) , po) = for any function /. 

Suppose then that po(y) / for some y £ M. Let us take f y (x) = 5(x,y)po(x) (i.e. we take 
the point source of the field po). Then f y (x) — (f) y = 5(x,y)p(x) — and 

(/«(*) " (fv) , Po) = P 2 o(v) - Poiv) (Po) = (19) 

which means po(y) = (po), and thus the equivalence class of p in Tq(C°°(M)/V) is zero, which 
makes the map injective. 

By the canonical transformation properties of 7^, one can extend this construction to arbi- 
trary (j). We have thus proven 

Proposition 1 The linear map given by 8 c %L{x) : Tq(C°°(M)/V) -»• Im(8 c %S) a C°°{M)/Nq 
where T$L(x) = T ( f ) S(x) — T ( f > (S(No))^fg, is a toplinear isomorphism for all (cj>,g,ir) provided 
ir ab £ 0. 

We have shown that it is a linear continuous bijection, and hence a topological linear isomor- 
phism. □. 

This singles out the equations we must solve for eft, namely, T^L = 0. But we can go 
further, not only can we form the Dirac bracket using {T,pL(x), n^y)} , but we can now 
use the implicit function theorem for Banach spaces for the function T^L^ _^ x ^ : T x C/V — > 

Im(5cT ( f > S) — C°°(M) /Nq. First we remind the reader of the formulation of the implicit function 
theorem for Banach spaces: 

Theorem 1 (Implicit function) Let X, Y, Z be Banach spaces. Let the mapping f : X xY — ?■ 
Z be continuously differentiable. If (xo,yo) € X x Y, such that f(xo,yo) = and d~Yf( Xo ,y ) '■ 
Y — > Z is a Banach space isomorphism from Y onto Z, then there exist neighborhoods U of xo 
and V of yo and a differentiable function g : U C X — > V CY such that f(x,g(x)) = and 
f(x,y) = if and only if y = g(x), for all (x,y) G U x V. 

In the above theorem, we substitute X H- T,Y H> C°°(M)/V,Z ^ C°°(M)/N , and / H> 
T(j,L. Now, for any (xo,yo) on the constraint surface {T^S = 0} n {D = 0} (which also implies 
that 7~y Q (L(xa)) = 0) we have that: 

Theorem 2 Around each point on the constraint surface there exist open subsets U x V £ 
r x C°°(M)/V such that there exists a unique (po : U — > C/V, such that the level surface formed 
from T<f,L = is given by 
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{(&y,7r <J ',0o[gy,7r y ],7r^) £ U x (T^v}- 

In other words, around each point on the constraint surface we can find the solution to Td,L(g, w, (ft, tt^) = 
for all (g, it, tt^) in a given open neighborhood by setting (ft = (fto® Furthermore, since we have 
in the process shown that dcT^L is surjective, the level surface formed by (T^L)" 1 (0) is regular. 
□. 

Thus we reach the surprising result that we can solve the second class constraints while 
eliminating the extra Stuckelberg variables! Although using the implicit function theorem in 
this context seems a bit contrived, we have a powerful theorem that is guaranteed to work 
whenever the operator A given in (7) is invertible. Although the operator A concerns only 
functionals in the original phase space V, the canonical transformation in the extended phase 
space not only tells us we can extend this result, but also that it will tell us a lot about the 
tangent of these maps on the (ft direction, enabling us in the end to use the implicit function 
theorem. 

We end this section by noting three things. First of all, the emerging theory, with constraints 
7~<j> (S (No)) , D , H is clearly not merely ADM in a CMC foliation, for it contains the extra, 
conformal constraint D, and the map 7^ , which do not appear in ADM under any guise. 
Secondly, we note that since No is not identically zero, the theory makes sense as a theory of 
geometrodynamics even if No(x) = for some x € S. This is why we have asserted that any 
ADM dynamical system has an SD counterpart, even if some vacuum ADM spacetimes do not 
have a CMC foiliation. Space-times not admitting a CMC foiliation are taken to be those where 
the lapse solution No reaches zero at some point, thus freezing time. In the geometrodynamical 
sense which we adopt here this assertion loses meaning. Lastly, we note that we did not solve 
the Lichnerowicz-York equation S(e^g, e 4< ^7r) = 0, where the </>'s are not restricted to maintain 
volume. In fact, we have not even solved our version T^S = 0, but only the non-homogeneous 
version T^L = allowing a global Hamiltonian to be left after reduction. 

In this direction we also note that one could possibly use the York-6'Murchada method of 
[6] to solve such types of equations, but the criteria for solvability becomes more complicated, 
and much harder to use. In [3], the authors originally set out to use both methods in parallel, 
but the flexibility and simplicity of the linear method in the end convinced us to leave the 
original York-O'Murchada method out. Furthermore, even if we were to use said methods, we 
must first come to the clean separation between first and second class constraints. I.e. we must 
come at least as far in our analysis as proposition 1 to know what is the equation we would like 
to apply a given method to, namely, what is the purely second class constraint in question. The 
advantage of our method is that from that point on it requires no extra work. 

3 Breaking the bound 

Now we mention briefly what becomes of our method once we couple different types of fields to 
the original construction. We will find that the equivalent resulting linear part of A is no longer 
positive, and we will find ways to circumvent this in the construction of SD. We have outlined 
the whole method of construction in the previous section so that we can follow a similar pattern 
for the more general case presented now. 

3.1 Coupling of matter fields 

The coupling of different types of matter fields was done more completely in [1], and reported 
more briefly in [7]. It was there found that if we assume an independent scaling of fields as in, for 
example in the scalar case, ift —> e^ift, for the propagation of gauge constraints develops 

6 Different solutions might exist around different open neighbourhoods, but we will not enquire on how these 
can be glued together. 
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pathologies which the authors could not amend. For this and other reasons, it was found that 
the natural choice for fields was "neutral coupling" . Neutral coupling is the particular choice of 
empty scaling a = 0, i.e. non-gravitational fields have conformal weight zero. This means that 
fields are only scaled in the sense that they are "carried along" by the scaling of the spatial 
metric. The method then proceeds as before, as we now briefly describe. 
We start with the general, first class constraints 

H(0 = Jd 3 x{7T ab (^g) ab + 7r A (L^ A )) 

S(N) = fd 3 x(^Tr ab G abcd ir cd -(R-2A)v^\+S m ^ (20) 
G a (X a ) = Jd 3 xG a (g ab ,^ A )X a ^g 

where we denote the "matter" (or rather, just non-gravitational) degrees of freedom collectively 
by i[) A and their canonically conjugate momenta tt a . We assume that the matter Hamiltonian 
H ma tter does neither contain ir ab nor any spatial derivatives of g ab . This is a reasonable assump- 
tion, and it holds for all fields known. We furthermore assume the constraint associated to 
internal gauge symmetries to be a functional of only the "position" variables ip A ,9ab- This is 
also a condition realized in all matter fields studied, and it simplifies treatment greatly. 

The Stuckelberg extension proceeds as before, but now we thus use the canonical transfor- 
mation generated by: 

F= [d 3 x (e^g ab n ab + ^ + (p A Ii A ) (21) 



which acts non-trivially only on the gravitational variables, as required by neutral coupling. 
The resulting canonical transformations leaves the extra constraint T&Q as in the vacuum case, 
the transformed diffeomorphism constraint H still weakly generates diffeomorphisms in the 
extended phase space (now also extended by the matter degrees of freedom), and the gauge 
constraint decouples from the conformal transformation, so that T$G a oc G a ~ 0. 
After imposing the gauge-fixing tt^(x) = we always get an equation of the form: 



{T (P S(N),tt4x)} = T (P 



--S(x) + 2^(A mMci N{x) - (A matter iV)) 



(22) 



If we here restrict to matter Hamiltonians which don't n contain spatial derivatives of the metric 
tensor nor metric momenta, we get: 

A — V/2 7r ^) p. 1 ( ^matter 3 \ 

^matter •— V ~ — ~ ti ~\ = g ab + - <3 ma ttcr \40) 

V5 2y/\g\ V °9ab 2 J 

where we note a slight deviation from previous notation, since here the subscript "matter" on 
the A operator serves to make it distinct from the vacuum case, and does not refer somehow 
only to the operator regarding the matter degrees of freedom. On the constraint surface T^S = 
and Q = 0, the end result is equivalent to taking 

A ~ 1 / \2\ °' afe(J afe 1 / ^matter 1 \ 

^matter « (V ~ — (Vf) j h — — — g ab - -i> matt er (24J 

12 g 2^/g V og ab 2 J 

Thus in this case the criterium for invertibility of the operator rests on: 

1 ( SS mBtteI 1 \ 1 2 . aab(T ab /r .t-\ 

2~fg {-Jg-r 9 ^ ~ 2 s —) * + — (25) 



the so called "bound" found in [3]. So the question that we face is: what happens if the bound 
is broken? 
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3.2 Finding the second class parts in the general case 

3.2.1 Properties of the solutions iVj 

Restricting S mMcr to not contain any metric derivatives, guarantees us that the symbol of the 
operator A mattcr given in (24) is defined by V 2 and thus A mattor is elliptic. Taken to be an 
operator between Sobolev spaces, A mattcr is then a Fredholm operator, and as such possesses 
the very important property that it has a finite dimensional kernel and cokernel [8], where 
the cokernel of a linear mapping of vector spaces / : X — > Y is the quotient space Y/im(f) 
of the codomain of / by the image of /. This immediately guarantees that such operators 
are invertible modulo compact operators. These are the technical facts that allows us to have 
a similar construction as in the invertible case. So let us follow the constructions from the 
previous section. 

First, to obtain the isomorphism we will still have to study the analogous equation to (10) 

A matter iV = ^A mattcr iv) (26) 

The difference is that besides one inhomogeneous solution AiV = c, we will have all the ho- 
mogeneous ones A mattcr iV = 0. Luckily, we know that A mattcr is Fredholm, and thus it has is 
a finite-dimensional kernel, Ker(A matter ). Let us assume that using a Gram-Schimdt algorithm 
(in the appropriate L2 norm) we can find an orthonormal set that generates Ker(A mattor ), we'll 
call them {iVj}j£/, where i £ N. Note that this is not the analogous condition required on the 
invertible case of the previous section. There we had to prove that it was possible to choose 
(Nq) = 1. Here, we are using (Ni,Nj) = 5ij, which requires no proof since we are assuming 
we are in a finite dimensional subspace of a Hilbert space. The reason for this departure will 
become clear soon. 

So any solution to the inhomogeneous equation has an added ambiguity given by a l Ni, 
where a 1 £ M, and assuming / has n elements, the space of solution of (26) is n + 1 dimensional. 
We note furthermore that since the solutions 2Vj depend on the point in phase space we are at, 
so does the dimension n. Lastly, we note that the different iVj's are linearly independent, and 
of course never identically zero as functionals of (g,ir): Ni(g,ir) ^ 0. 

3.2.2 Separating out the purely second class constraints 

Again, we have (11) 

TS(x) : r x T*(C°°(M)/V) x T mattcr -> C°°(M), (27) 
where r mattcr denotes the phase space of the matter fields. Thus, following (14): 

ScTZSbJo, Ko,ir%^,A)\<t>=o ■ T (C°°(M)/V) -> C°°(M) (28) 
and finally from (16) 

{5 C %SY ■ N = A matter iV - (A matter iV) (29) 
and from the eliipticity we in fact have 

n 

hn(5 e %S) N t = C°° (M) (30) 
Now we construct the tentative operator 

n 

L (x):=S(x)-^g-J2 S ( N ^ ( 31 ) 

i=0 
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We can readily check that the tangent map Sc%f,Lo^ _ indeed stays (weakly) the same, but 

{Lo,a l Ni) 7^ 0, even if we were to assume the conditions (iVj) = 1. At this point it is already 
clear that it will be hard to recover the original version of Shape Dynamics, as indeed we don't. 

Thus, we are forced to depart the usual construction of Shape Dynamics, and now we use 
the fact that we have made the set of generators iVj orthonormal. The operator we are to use 
becomes: 

n 

L(x) := S(x) -^Y, S ( N i) N i( x ) (32) 

i=0 

Now it becomes clear that (L,Nj) = 0, and the tangent map o~cTtf>L^ =0 is still (weakly) the 
original. Thus 

(T^L)^ =f{x) : r x C°°(M)/V lm(S c T^S). (33) 

To finish the proof that 8cTd,L is a topological linear isomorphism when restricted to the 
correct space, we still need to check whether the space of conformal factors needs restrictions 
so that the resulting operator is injective. A priori, we would expect that we would also have 
to quotient the space of conformal factors by an n-dimensional space. This would mean that 
there would be an (n + 1) x (n + 1) set of generating Hamiltonians. As we will see, surprisingly 
this is not the case. 

Proposition 2 

(ImiScT^S)*) 1 - = Ker(5 c %S) = V 
where V is the space of constant functions. 

The proof we offer is not very clean, and there might be better ways of attaining the same result. 
In any case, the proposition yields the surprising result that we can still completely solve the 
second class constraints for the conformal factor <j) , as in the vacuum case. 

Again we must analyze (ImdcT^S ) . To provide some guidance, and also a necessary pre- 
liminary, 7 consider the following simplified case: suppose that (ScT^S)* was of the simpler form 
(ScT^S)* • N = A matter iV\ Of course, A matter is self-adjoint, so we know that Ker^cT^S)*) = 
Ker(5cT^S) and we already know the complete characterization of the decomposition once we 
know the kernel. Nonetheless, let us follow through the computations. Using the L2 inner prod- 
uct structure we can thus easily parametrize the image of (5cT<pS)* which is just (Ker(5cT < i } S)) ± , 
as: 

{/(*) " E </' N *> I / G C°°(M)} (34) 

i 

Thus we have that p G Kev{5 c %S) if and only if, for all / G C°°(M): 

( P J-^2(f,N l )N i ) = (35) 

i 

By choosing, in the same way as done in (19), for a given point y G M (for which p{y) 7^ 0), 
f y (x) := p(x)5(x,y), we get: 

p(y)-J2N i (y)(p,N i )=0 (36) 

i 

which is a consistent equation solved by p(x) = Yli at Ni(x), as expected. 

Now, let us move forward to the full problem. This time we are trying to find the space 
orthogonal to 

{A matter iV(x) - (A matter iV) I N G C°°(M)} (37) 

7 This is the second warm-up to the problem. The first one was given in section 2.2.2, see equation (19). 
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We know that we can further parametrize this by substituting 

A matter iV(x) -> f(x) - E (f, Ni) Ni(x) 

i 

using self-adjointness and knowledge of the kernel of A matter . We thus have that the condition 
for p to be in the kernel of ScT^S is that, for all / £ C°°(M): 

</, p) - E (/- N i) (p> N i) - (/) (p) ^ + E </> ^> w> (p) ^ = ( 38 ) 

Contrary to what we might have supposed, p = a l iVj is not a solution for whichever set of 
ctj's we choose, since then (38) is equivalent to 

(f)-Y / (f,N i )(N i ) = (39) 

i 

and {Ni} only forms a basis of a finite-dimensional subspace (thus Y2i •) Ni is not a decom- 
position of the identity). We can also easily see that for k = const, p{x) = k is a solution of 
(38), and thus for p(x) = h(x) + k, h(x) has to separately be a solution. 

Finally, choose in (38) f(x) = p(x)5(x,y), as before. After a tiny bit of algebra we get: 

p(x) = (p) + Ni(x)({p, Ni) - (Ni) (p) V) (40) 

Of course, this implies that p(x) = h(x) + k, where h(x) = ^ a J iVj(x) for some choice of a* 6 R. 
As we saw, h{x) must be separately a solution, but from (39), we know that it isn't. Therefore 
the only solution is a 1 = 0, thus p(x) = k . □ 

This allows us to find the first and (purely) second class constraints, a necessary step for 
our construction. From here we can know what is the equation we have to explicitly solve 
(i.e. the purely second class constraints). Hence even if we are to use the York-O'Murchadha, 
Leray-Schauder method, we must come this far to know what is the equation we would like to 
apply it to. That said, we have been unable to work out how said method might work for the 
second class constraint found here. 

By the canonical transformation properties of 7~d>, one can extend this construction to arbi- 
trary <j). We now have 

Proposition 3 At the constraint surface, the linear map given by 

n 

5 c %L(x) : T (C°°(M)/V) -> Im{5 c %S) ~ C°°(M)/ 0iV, 

8=0 

where 

n 

T ( pL(x) = T<pS(x) - E TtiSiNMixfiy/g (41) 

i=0 

is a toplinear isomorphism. 

We have shown that it is a linear continuous bijection, and hence a topological linear isomor- 
phism. □. 

3.3 The reduced system 

Proposition 3 tells us what part of the constraint is "purely second class" and thus what part 
of the constraint should be explicitly solved for. Again we have the Dirac bracket and use the 
implicit function theorem for Banach spaces for the function T^L around the constraint surface 
to assert that 
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Theorem 3 There exists a unique <fio : U C (T x T mattcr ) — > V C (C°°(Af)/V), suc/i t/iai the 
level surface formed from T^L = is given by 

{(gij,TT ij \vA^A,0obij,^ J '],^) £ u x (r )| y }. 

In other words, given any point on the constraint surface, the solution to T c f,L(g,TT,(j),-K ( f ) ) = 
for all (g,TT,ip A ,7TA,^(j,) on a given open set around this point is obtained by setting cf> = <pQ. 
However, even though we still have that for each point b~cT<hL is surjective, the level surface 
formed by (7^ ) L)- 1 (0) may not be regular, since the number of linearly independent solutions 
Ni may change from point to point in phase space. This fact will not be dealt with here. 

We also note that one of the initial fears of breaking the bound, namely that there would be 
a different 4> for each homogeneous solution, is not realized, since (ft comes from the solution 
of T<f>L = 0, and not from the "lapse fixing equation". The local implication that there might 
be different <fi functionals around different points (g, ir, tp A , tta) also does not concern us here. 

Finally, we deal with the reduced system. First, note that for none of the solutions Ni, we 
have that iV, = 0, and the N are furthermore linearly independent. This entails that if we 
have a constraint of the form a % Ni = 0, we can only satisfy it (remember the constraints 
must be satisfied for all x € M) by setting a 1 = 0. Combining this with the result that 
T^f(g,TT,TT A ,ipA)(x) = if and only if f(g, n, ir A , ip^)(x) = (see the appendix in [2]), we have 
that the reduced constraints 

n 

^T^iSmNiix)) (42) 

i=0 

are equivalent to the n + 1 independent constraints: 

{T^SiN)) \i = 0,...n} (43) 

This resolves a possible paradox with the set of leftover physical degrees of freedom of the 
theory, which would have been enlarged if there was no such simplification of the constraints. 
We note that with this new method, which reduces to the vacuum case, previous considerations 
about iVo (such as the possibility of choosing it such that (Nq) = 1) are turned obsolete. 
Finally, we arrive at the generalized Shape Dynamics total Hamiltonian: 

#sd = a l T^(S(Ni)) + J d 3 x (> 6 (L^U + tt a (C^) a + e k ^G a (g ab , MKVd) ■ (44) 

where we have used the decoupling of the Gauss constraints related to the extra fields, T^G a = 
e m<t> G * which 

appears for some m for all fields studied (for example, for Yang- Mills m = 6). 

We take the opportunity here to note the important fact that this theory is not the same 
as the original SD, in the sense that it does not reduce to SD even when there are no non- 
trivial homogeneous solutions to A mattcr . What is different is the form of T$L(x) given in (41). 
This does not reduce to T^S{x) — 7^,(5 '(A^)) ^fg , the original case realized for matter obeying 
the bound. Thus, although the form of the global Hamiltonians is the same, {T,f )o (S(Ni)), the 
functional (fr might be different, even if the bound is satisfied. 

In any case, it is a trivial exercise to check that the non-zero part of the first-class constraint 
algebra is given by: 

{H a (r) a ),H b (Z b )} = H a {[lfJ\ a ) 

{H a (C a ),D(p))} = D{C^p) (45) 
{a% o (S(Ni)),P% (S(Ni))} = a i FH a (N j V a Ni-N i V a N j ) (46) 

(47) 
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4 Conclusion 



The present work considers the question put forward previously in work coupling different fields 
to Shape Dynamics [3, 7]. Namely, it was found that for some types of couplings, including the 
mere inclusion of a cosmological constant, certain bounds on the previous construction of SD 
were imposed. In usual GR language, these bounds ensured that there existed a unique lapse 
propagating a CMC condition, and were technically implied by the triviality of the kernel of a 
given differential operator, called here A matter . Thus not only for completeness of the theory, 
but also if we would like to make statements about cosmology in SD, it was seen as important 
that we sorted out such limitations in the theory, which is what this paper was aimed at. 

Using the same method of "excising" kernels and co-kernels to obtain isomorphisms as 
brackets of "purely second class constraints", a method fully described in section 4.3 of [2], 
we attempted to extend previous results of standard SD beyond such bounds. Surprisingly, 
we found a theory slightly different than SD, even when the bounds are respected. In such 
case the distinction will be manifested only by the appearance of a different scalar functional 
(j) [g,7r,ip A ,tta]{x), appearing in a global Hamiltonian of the same form as the original. The 
other constraints and the constraint algebra are unchanged. The general theory found here still 
exists as a conformally invariant theory, possessing all the advantages of the original version, but 
now it may have also a finite number of weakly commuting global Hamiltonians. The number 
of Hamiltonians is given by n + 1, where n is the dimension of the kernel of A matter , and yields 
the number of different lapses propagating a CMC condition. 

These results also serve to make more distinct the role developed by cj) and the "lapse 
solutions" Ni, which were previously thought to be 1-1, or correspondent. Before these results 
we believed that for each homogeneous solution of A mattcr there would also exist a corresponding 
different solution of the second class constraints for the Stuckelberg field <fi. This would yield a 
space of (n + 1) x (n + 1) additional constraints which might have turned out to be inconsistent. 
Had there been different reductions, i.e. different solutions for </>, interpretation of the result 
would have been extremely difficult, since we would have had several generators of evolution 
which could not have been shown to commute. To be more explicit, departing from {T^f, Tsh} = 
T,j){f, h}, for a single reduction the commutator will be simple to compute, but if we had two 
reductions 4>i,4>2 it becomes unclear what the result might have been. However, the reduction 
process was shown to be unique, and these fears were not realized. This surprising result 
was attained in proposition 2. This work gives us the picture of a general theory of Shape 
Dynamics, always with a finite number of generators of time evolution, with a related freedom 
to continuously shuffle between them. 

We would also like to highlight the use of the implicit function theorem as a purely technical 
tool. Coupled to canonical transformations, it is able to yield local solvability of equations of 
Lichnerowicz-York type from Poisson brackets in the original phase space T x r mattcr . This is 
a very different approach than Leray-Schauder theory of non-linear differential equations poly- 
nomial in (f>, employed by York-O'Murchadha [6] to solve a "conformalized" scalar constraint. 
The present method yields a simple, linear criterion for local solvability of such equations. As 
far as the author knows, this technique has not been previously employed. More importantly, 
to find the equations one needs to solve for, one must separate out the "purely" second class 
constraints. I.e. those that we would like to solve as definitions of the auxiliary variables. This of 
course has already happened in the vacuum case, whence we arrived at the "non-homogeneous" 
version of the LY equation: k[g, ir] = T^S^x), where k[g, n] was a spatial constant. It is only 
from this point on that one can even discuss using the York-O'Murchadha method. Thus even 
to formulate the equation one wants to solve in terms of the 4> field (e.g.: k[g, ir] = T<pS(x) ) 
one must determine a split of the original scalar constraint such that one of the parts has in- 
vertible Poisson brackets with the CMC constraint. This is a necessary technical point that our 
dynamical study must arrive at in any case, and from here on the use of the implicit function 
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theorem requires no extra work. 

In connection to this, let us mention a practical difference between the generic case presented 
here and the vacuum case. In the vacuum case, the fact that we could find the solution <\> by 
use of the equation k[g, tt] = %pS(x), where k[g,ir] was a spatial constant, was due to the 
fact that the purely second class part of T^S was of the form T c j ) S{x) — T<pS(No). Now the 
correspondent second class part, given by equation (41), is not nearly as simple, and the rhs is 
spatially dependent. In this case the York-O'Murchadha method does not yield an alternative 
method of proof, it is only through the implicit function theorem that we can find the right 
theory. 

Regarding actually finding solutions, in the general case then one has no choice but to find 
the individual solutions N{, plug them in and work out solutions of <j) . Of course, it remains 
true that if (g, n) belong to a physical ADM solution S(g, tt) = 0, then we still obtain <j) = 0. 
Thus a useful approach to find non-trivial solutions of SD is, as in the vacuum case, to perform 
a perturbative analysis around a given ADM spacetime. 
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